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A continuum extrapolation of static four- and two-quark energies calculated in quenched SU(2) is done based on 
Sommer's method of setting the scale. A model for four-quark energies with explicit gluonic degrees of freedom 
removed is fitted to these energies and the behavior of the parameters of the model is investigated. 



1. Introduction 

Our group is attempting to understand hadronic 
interactions from first principles by developing 
a potential model for multi-quark systems with 
explicit gluonic degrees of freedom removed (0] 
and references therein). Such a model would 
shed light on multi-quark bound states and e.g. 
meson-meson scattering. 

We have been simulating systems of four static 
quarks in quenched SU(2) at (3 = 2.4 with a 
16 3 x 32 lattice and the Wilson action. Various 
geometries such as rectangular, linear and tetra- 
hedral have been considered to get a set repre- 
senting the general case. 

To find the lattice artifacts in energies and their 
parameterization, I have performed continuum 
extrapolations based on Sommer's method of set- 
ting the scale [§ at (3 = 2.35, 2.4 (16 3 x 32), 2.45 
(20 3 x 32), 2.5 (24 3 x 32) and 2.55 (26 3 x 32) for 
squares and tilted rectangles, since these geome- 
tries exhibit the largest binding energies. 



2. A Model for Four-Quark Energies 

The model is based on static two-quark poten- 
tials for different pairings of the four-quarks. In 
the case of two pairings or basis states A,B, the 
eigenenergies \ are obtained by diagonalising 



with 
N = 

V = 




(2) 



[V - AjN] ^ = 0, 
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(1) 



where Vij represents the static two-body poten- 
tial between quarks i and j and Vab is from the 
perturbative expression Vij = —Af(N c )Ti ■ TjWy. 
The normalization is chosen to give for a color 
singlet state < [ij] \Vij\[ij] >= v tj . The 
four-quark binding energies Ei are obtained by 
subtracting the internal energy of the basis state 
with the lowest energy, e.g. Ei = A-j — (fi3 + ^24)- 

A central element in the model is the phenomeno- 
logical factor / appearing in the overlap of the 
basis states < A\B >= f/N c for SU(N C ). 

• / is a function of the coordinates of all 
quarks, making the j^-Vab in eq. || a 
four-body potential 

• in the weak coupling limit / = 1 

• a parameterization / = f ce ~k A b s A-k P VbsP ; 
where A, P are the minimal area and 
perimeter, 65 the string tension and 
f c ,ltA, kp are to be fitted 

Perturbative calculation to 0(a 2 ) produces the 
two-state model of eq. with / = 1 || . 

The above model using only the ground state of 
the two-body potential works for rectangles but 
fails for tetrahedrons. 
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3. Setting the Scale 

Sommer || has designed a popular new way to set 
the scale that avoids the long-distance limit in the 
definition of the string tension 65 = linv^oo F(r). 
After calculating the force F(r/a) between two 
static quarks, the equation 

(r /a) 2 F(r /a) = c (3) 

is solved with c = 1.65 for ro/a. Equating to the 
continuum value ro ~ 0.5 fm gives a. 

Choosing c = 1.65 produces agreement on ro = 
0.49 fm for the Cornell [0] and Richardson || non- 
relativistic effective potentials, while other mod- 
els §,0] give ro from 0.44 to 0.56 fm. However, 
setting c = 2.44 makes three potential models 
agree on ro = 0.66 fm while others give 0.625 or 
0.64 fm. 

Runs were performed for square and tilted rectan- 
gle (TR) geometries, with resulting scales shown 
in table 0. Here a is calculated using c = 1.65, 
a 11 using c = 2.44 and a bs from the string ten- 
sion using a continuum value y/b~s — 0.44 GeV. 
We can see that at each j3, a bs /a w 1.09, and be- 
comes « 1 if we set V&s = 0.478(4) GeV. Mean- 
while a 11 needs no such adjustment to agree with 

So it seems c = 2.44 is a better choice than 
c = 1.65. I chose the former value, keeping in 
mind that the difference can be accidental due 
to the significant uncertainties in the nonrela- 
tivistic potential models and the difference of 
the quenched SU(2) string tension from a phe- 
nomenological value. 

4. Extrapolating Two-Body Potentials... 

To extrapolate, we need values of energies at dif- 
ferent P's corresponding to the same physical size. 
Sample two-body extrapolations (from a total of 
29 potentials) are shown in fig. for potentials 
with separation from 0.17 to 0.52 fm involved 
in squares. Other potentials look similar. A 
quadratic fit is clearly preferred. 



p 


a [fm] 


n — ft — t 

a 11 [fm] 


a" s [fm] 


2.35 


0.1302(6) 


0.1401(14) 


0.1408(11) 


2.35(TR) 


0.1299(7) 


0.1397(6) 


0.1402(5) 


2.4 


0.1101(10) 


0.1186(10) 


0.1194(9) 


2.4(TR) 


0.1103(5) 


0.1190(4) 


0.1200(11) 


2.45 


0.0918(7) 


0.0997(8) 


0.0997(7) 


2.45(TR) 


0.0929(7) 


0.1002(7) 


0.1010(6) 


2.5 


0.0793(5) 


0.0856(9) 


0.0866(9) 


2.5(TR) 


0.0789(5) 


0.0851(6) 


0.0859(6) 


2.5(W) 


0.0793(3) 


0.0856(6) 


0.0864(5) 


2.55 


0.0675(5) 


0.0730(5) 


0.0741(7) 


2.55(TR) 


0.0673(9) 


0.0727(10) 


0.0734(8) 



Table 1 



Values of a for each j3. (W) is for Wuppertal. 



5. ...and Four-Body Binding Energies 

Binding energies for the simulated geometries can 
be interpolated using the /-model. I used the pa- 
rameterization / = f c e~ kAbsA , as the perimeter 
term was not needed. Results for squares with the 
length of a side from 0.17 to 0.52 fm are shown 
in fig. 0. Results for tilted rectangles look sim- 
ilar. The weaker dependence on (ro/a) 2 shows 
that there are significantly less lattice artifacts 
than for the two-body potentials. 

6. Extrapolating the Parameters of the 
Model 

Values of the parameters were obtained by fit- 
ting / = f c e~ kAbsA to the energies of squares and 
tilted rectangles together at each (3. The normal- 
ization f c and the constant multiplying the area, 
kA, are extrapolated in fig. 0. 

The parameters are changing quite a lot, while / 
itself is more constant at different (3's. This can 
be taken as a sign of a poor parameterization. 

7. Continuum Fit of the Model 

After extrapolating the two- and four-body ener- 
gies the /-model was fitted to the continuum val- 
ues. Both linearly and quadratically extrapolated 
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0.020 „ 0.030 

(a/r ) 2 




0.010 0.020 ,0.030 0.040 0.050 

(a/r ) 2 



Figure 1. Two-body potentials for squares with 
length from 0.17 (top) to 0.52 fm (bottom) 



Figure 2. Four-body binding energies for squares 



energies were fitted. Results of the most relevant 
fits and parameter extrapolations are shown in 
table § 



type 


fc 


k A 


x 2 


lin/quad 


1.04(4) 


1.04(6) 


0.7 


lin/quad, f c = 1 


1 


0.98(3) 


0.8 


quadratic 


1.08(8) 


1.37(12) 


0.4 


quadratic, f c = 1 


1 


1.26(6) 


0.5 


linear extrap. 


1.04(3) 


0.83(6) 


2.6/3.3 


quadratic extrap. 


1.10(5) 


1.09(11) 


2.7/0.4 



Table 2 

Continuum fit and parameter extrapolation re- 
sults. Fit type "lin/quad" denotes linear extrap- 
olations for four-body energies and quadratic for 
two-body potentials. 



Our best estimate for k A is slightly above one, 
while f c can be safely set to one. 

In the strong coupling limit a factor exp (—bsA), 
similar to /, appears in the off-diagonal elements 
of the Wilson loop matrix, A being the minimal 



surface bounded by straight lines connecting the 
quarks. It has been argued that weaker couplings 
lead to a smaller transition area J8|. If this is 
the case, it is not reflected in our best estimate 
of k A - The larger transition area pointed to by 
kA > 1 could be explained by the finite width of 
flux tubes as compared to the lines in the strong 
coupling approximation. 

Continuum fits were also performed for squares 
extrapolated using the original value c = 1.65 in 
equation ||. The resulting values of f c and fc^ for 
fits to all linear and quadratic data were within 
errors with fits to extrapolations using c = 2.44. 

8. How to develop the model? 

The failure of the simplest form of the /-model 
to predict the binding energies of the tetrahedral 
geometry has been proposed to be due to a depen- 
dence of the four-body ground state energy E± on 
the excited state of the two-body potential. 
This was investigated using three fuzzing levels, 
all with c = 4, producing different ground- and 
excited state overlaps of the two-body paths. 

• fuzz level 150: overlaps 97.9-98.7 %, inter- 
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Figure 3. f c and fc^ continuum extrapolations 



mediate convergence of £4 

• fuzz level 20: overlaps 99.1-100 %, best con- 
vergence of £4 

• fuzz level 0: overlaps 85.3-98.2 %, worst 
convergence of £4 

Thus the excited state contribution to £4 comes 
from the overlap of the ground state of a gluon 
field between two quarks with the excited state of 
a gluon field of another quark pair. In ref. j^] this 
is investigated further and a generalization of the 
model is introduced that is capable of reproducing 
the degenerate ground state energy of a regular 
tetrahedron. 

It is evident that using a high fuzzing level 150 
(as in Wuppertal [|o|) does not necessarily lead 
to better ground state overlaps on a 24 3 spatial 
lattice. 

The effects of instantons are another candidate 
for a missing feature in our model. Instanton 
liquid models include neither perturbative one- 
gluon exchange (Coulomb part of the two-quark 
potential) nor confinement (linear part of said 
potential). They still predict the nucleon to be 
bound with the right mass O]. Our model gets 



the energy scale only from the two-body potential 
and the overlap factor / can hardly account for 
instanton effects. 

Our next project is an attempt to uncover the 
microscopic origins of the simulation results by 
relating the model to flux distributions via energy 
sum-rules similar to those derived by C. Michael 
@- 
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